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Abstract 


Given  the  shortest  Hamiltonian  path  (or  tour)  in  an  undirected 
weighted  graph,  the  sensitivity  analysis  problem  consists  in  finding  by  how 
much  we  can  perturb  each  edge  weight  individually  without  changing  the 
optimality  of  H°. 

The  maximum  increment  and  decrement  of  the  edge  weight  that  preserve  the 
optimality  of  H°  is  called  edge  tolerance  with  respect  to  the  solution  H°.  A 
method  of  computing  lower  bounds  of  edge  tolerances  based  on  solving  the 
sensitivity  analysis  problem  for  appropriate  relaxations  of  the  shortest 
Hamiltonian  path  (tour)  problem  is  presented. 
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1.  Introduction 

^  The  problem  considered  in  this  paper  belongs  to  so  called  sensitivity 

r 

analysis  in  combinatorial  optimization  (see  e.g.  [2]).  This  term  is  used  for  a 
phase  of  solution  procedure  when  an  optimal  solution  of  problem  has  been 
already  found  and  additional  calculations  are  performed  in  order  to 
investigate,  how  this  optimal  solution  depends  on  changes  of  problems 
parameters.  ^ 

In  this  paper  two  well  known  (see  e.g.  [7])  combinatorial  optimization 
problems  are  considered:  the  shortest  Hamiltonian  path  problem  in  undirected 
weighted  graph  and  the  symmetric  traveling  salesman  problem.  It  is  assumed 
that  an  optimal  solution  of  given  problem  is  known.  The  goal  of  sensitivity 
analysis  consists  in  finding  by  how  much  we  can  perturb  each  edge  weight 
individually  without  changing  the  optimality  of  the  solution.  The  maximum 
increment  and  decrement  of  the  edge  weight  that  preserve  the  optimality  of 
solution  are  called  the  edge  tolerances  with  respect  to  this  solution. 

In  this  paper,  a  method  of  computing  lower  bounds  of  the  edge  tolerances 
with  respect  to  the  optimal  solution  of  the  shortest  Hamiltonian  path  problem 
and  traveling  salesman  problem  is  described.  The  method  is  based  on  solving 
the  sensitivity  analysis  problem  for  appropriate  relaxation  of  the  original 
optimization  problem.  >  A  general  idea  of  this  approach  was  presented  in  [8]. 
In  this  paper  we  give-^a  description  of  the  approach  and  its  microcomputer 

l 

implementation  and  we  report  preliminary  results  of  computational  experiments. 

This  paper  is  organized  as  follows.  In  section  2  we  introduce  a  notation 
and  give  some  preliminary  results  concerning  the  relations  between  the 
sensitivity  analysis  for  the  original  problem  and  its  relaxation.  In  section 
3  we  describe  algorithms  for  performing  a  sensitivity  analysis  for  problem 
relaxations.  A  choice  of  appropriate  relaxation  of  the  original  problem  is 


discussed  in  sections  3  and  4.  Section  4  contains  also  a  description  of 
implementation  of  the  method  and  results  of  numerical  experiments. 


2.  Notation  and  preliminary  results 

Let  G  =  (V,E,C)  be  an  undirected  weighted  graph  with  a  set  of  vertices 


V  =  {1,...,n}  and  a  set  of  edges  E  =  {e^,...,e^}  c  V  x  V. 


C  «  R 


where  R  =  R  u  {»} ,  is  a  matrix  of  edge  weights.  (If  e  =  (i,j)  ^  E,  then 

c(i»J)  =  “•)  The  subgraph  (V,Q,C)  of  G  will  be  identified  with  a  set  of  its 

edges  Q  and  by  1(Q)  =  J  c(e)  we  will  denote  a  weight  of  the  subgraph. 

esQ 

Let  H  be  the  set  of  Hamiltonian  paths  in  G  with  fixed  ends  in  vertices 
1,  n  and  let  H  denote  the  set  of  Hamiltonian  tours  in  G.  Two  well  known 
combinatorial  problems:  the  shortest  Hamiltonian  path  problem  (SHPP)  and  the 
traveling  salesman  problem  (SHIP)  are  formulated  as  follows 


min{l(H)  :  H  «  H} 


min{l(H)  :  H  «  H} 


(SHPP) 


(SHIP) 


In  this  paper  the  shortest  Hamiltonian  path  problem  will  be  mainly 
considered.  The  approach  for  the  traveling  salesman  problem  is  similar;  the 
differences  are  pointed  out  if  necessary. 

Assume  that  H°  is  a  (known)  optimal  solution  of  the  SHPP  in  the  graph  G, 

i.e. , 


H°  =  arg  min{l(H)  :  H  «  H} 


The  tolerance  problem  is  formulated  as  follows: 
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Given  find  for  e  «  E  values  c'*"(e),  c~(e),  such  that  is  optimal 
for  any  perturbed  graph  G'  =  (V,E,C’),  in  which  c'(i,j)  =  c(i,J)  if 
(i,J)  i  e  and  c(e)  -  c“(e)  <  c*(e)  <  c(e)  +  c'^(e). 

The  values  c'''(e),  c~(e)  are  called  upper  and  lower  tolerances  of  the  edge  e  with 
respect  to  the  optimal  solution  H°.  Edge  tolerances  with  respect  to  optimal 
solution  of  the  SHIP  are  defined  in  the  same  way. 


Let 


H  =  {H  «  H  :  e  «  H} 
e 


and 


f/®  =  {H  «  H  :  e  H} 

The  following  proposition  expresses  the  edge  tolerances  c''’(e),  c~(e),  e  s  E, 
by  auxiliary  optimization  problems  over  sets  H®,  (We  will  assume  that  if 

a  minimization  problem  is  infeasible,  then  its  optimal  value  is  equal  to  ”) . 


Proposition  1.  If  e  «  then  c~(e)  =  “  and 

c'*’(e)  =  min{l(H)  :  H  s  H®}  -  1(H°). 

If  e  ^  H°,  then  c‘^(e)  =  ®  and 

c~(e)  =  min{i(H)  :  H  «  ff  }  -  1(H°). 

e 


(1) 


(2) 


a 


I 


I 
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Proof.  Consider  an  edge  e  «  H°.  It  is  obvious  that  any  decrement  of  the 
weight  c(e)  does  not  change  the  optimality  of  H°,  so  c~(e)  =  ®.  If  the 
weight  of  e  increases  and  the  weights  of  all  other  edges  remain  unchanged, 
then  weights  of  all  Hamiltonian  paths  belonging  to  Hg  also  increase  in  the 
same  way,  but  weights  of  paths  in  H®  are  still  the  same.  Therefore  H° 
remains  optimal  as  long  as  the  increase  of  the  weight  of  e  is  not  greater  than 
the  difference  between  the  weight  of  the  shortest  Hamiltonian  path  in  H®  and 
the  value  1{H°).  The  proof  of  the  second  part  of  Proposition  1  is  analogous. 

D 

Similar  fact  may  be  proved  for  edge  tolerances  in  the  SHTP. 

Proposition  1  suggests  that  a  calculation  of  edge  tolerances  may  be  a 
difficult  task,  because  in  order  to  find  the  tolerances  for  a  particular  edge, 
one  has  to  know  the  optimal  value  of  an  auxiliary  optimization  problem,  which 
is  in  general  as  difficult  as  the  original  SHPP  (unless  this  value  is  a  by¬ 
product  of  solving  the  original  problem).  Another  explanation  of  difficulty 
of  this  sensitivity  analysis  arises  from  the  observation  that  the  tolerance 
problem  is  closely  connected  to  a  problem  of  finding  adjacent  vertices  in  the 
SHPP  or  the  SHTP  poly tope,  which  is  known  to  be  NP-hard  [7]. 

The  goal  of  this  paper  is  to  propose  an  approach  which  allows  to  compute 
in  an  efficient  way  lower  bounds  of  edge  tolerances,  i.e.,  values  d'^Ce), 
d“(e),  e  s  E,  satisfying  the  conditions  d''’(e)  c'*’(e),  d“(e)  <  c"(e), 

e  «  E.  Such  lower  bounds  are  also  of  practical  value,  because  they  imply 
that  for  particular  edge  e,  the  solution  H°  remains  still  optimal  if  the 
weight  of  e  belongs  to  an  interval  [c(e)  -  d”(e),  c(e)  +  d‘*‘(e)].  Calculation 
of  lower  bounds  seems  to  be  much  easier  than  calculation  of  edge  tolerances, 
because  in  order  to  find  c'*'(e),  c~(e)  one  must,  in  fact,  exploit  necessary  and 


sufficient  conditions  of  the  optimality  of  H°.  To  calculate  d‘^(e),  d“(e)  it 
is  enough  to  have  only  some  sufficient  optimality  conditions  for  H°.  The 
notion  of  necessary  and  sufficient  optimality  conditions  in  seldom  the  case  in 
combinatorial  optimization,  whereas  sufficient  optimality  conditions  are 
provided  by  different  relaxations  of  the  original  problem  and  related  dual 
problems.  The  choice  of  appropriate  relaxation  is  discussed  in  sections  4  and 
5.  In  this  paper  as  a  relaxation  of  the  SHPP,  the  shortest  spanning  tree 
problem  (SSTP)  is  chosen,  and  to  calculate  bounds  of  edge  tolerances  for  the 
SMTP,  the  shortest  1-tree  problem  (S1TP)  is  used  (see  e.g.  [7]). 

Let  us  consider  a  pair  of  problems  -  the  SHPP  and  the  SSTP  -  and  let 
v(SHPP),  v(SSTP)  denote  its  optimal  values,  i.e., 

v(SHPP)  =  rain{l(H)  ;  H  «  //} , 

v(SSTP)  =  min{l(T)  :  T  «  T}, 

where  T  is  a  set  of  spanning  trees  in  G.  Usually,  the  SSTP  is  not  a  good 
relaxation  of  the  SHPP  (if  we  measure  a  quality  of  relaxation  by  the 
difference  between  the  optimal  values  of  both  problems).  But  it  is  well  known 
that  this  difference  may  be  significantly  reduced  (see  e.g.  [7],  Chapter  10) 
by  appropriate  modification  of  edge  weights.  This  modification  consists  in 
replacement  of  the  original  edge  weights  c( i , j ) ,  ( i ,  J )  s  E,  by  values 
cP(i,j)  defined  as  follows: 


where  p(i),  p(J),  i»J  «  V,  are  elements  of  so  called  penalty  vector 

p  =  (p( 1 ) , • . . ,p(n) ) ^  «  r".  Denote  by  cP  modified  edge  weight  matrix  and  let 
GP  =  (V,E,cP).  The  weight  of  subgraph  Q  in  qP  will  be  denoted  by  lP(Q).  It 
is  well  known  that  the  such  modification  of  the  graph  does  not  change  the  set 
of  optimal  solutions  of  the  SHPP.  The  following  proposition  states  that  this 
is  also  true  for  edge  tolerances.  The  same  facts  hold  also  for  the  SMTP. 

Proposition  2,  Edge  tolerances  c'*’(e),  c“(e),  e  «  E,  are  the  same  for  any 
modified  graph  gP  =  (V,E,cP),  p  «  r’’. 

Proof.  This  is  a  simple  consequence  of  Proposition  1.  It  is  easy  to  see  that 
for  p  fi  R*^  the  value  d(H',H")  =  iP(H')  -  lP(H'')  does  not  depend  on  p  for  any 
H',  H”  «  H.  But  according  to  (1)  and  (2),  if  c'*’(e),  c~(e)  <  »,  then 
c'''(e)  =  d(H®,H°),  e  «  E,  ancj  c"(e)  =  d(H-,H°),  e  «  E\H°,  where 
H®  =  arg  min  {1(H)  :  H  «  H®},  =  arg  min{l(H)  :  H  s  H^} .  If  for  some 

e  «  c^(e)  =  ®  or  for  e  «  E\H°,  c~(e)  =  this  means  that  corresponding 

set  H®  or  is  empty;  which,  obviously  does  not  depend  on  the  vector  p. 

D 

Let  p  «  r”  be  arbitrary  penalty  vector  and  define  A(p)  to  be  equal  to 
the  difference  between  the  optimal  values  of  the  SHPP  and  the  SSTP  in  gP. 
Moreover,  let  tP  be  the  optimal  solution  of  the  SSTP  for  gP  and  define 

tp(e,TP)  (tp(e,T*^)),  e  «  E,  be  an  upper  (lower)  tolerance  of  e  with  respect 

to  tP  regarded  as  an  optimal  solution  of  the  SSTP  in  gP,  i.e., 

tp(e,T*^)  (tp(e,TP))  is  equal  to  the  maximum  increment  (decrement)  of  the 

weight  of  e,  which  does  not  change  the  optimality  of  tP.  Then  the  following 


fact  hold: 


Lemma  1:  For  p  s  R"  and  e  «  H“nT^u  (E\H^)  n  (E\T»') 


c‘*’(e)  >  tp(e,TP)  -  A(p) 


c'(e)  >  t"{e,TP)  -  A(p) 


Proof.  We  will  prove  only  (4);  the  proof  of  (5)  is  analogous.  If  e  s  E\H  , 
then  c^(e)  =  ®  and  (4)  holds.  Assume  then  that  e  e  H°riT^  and  let 


t®  =  min{lP(T)  :  T  s  T®} 
P 


where  =  {T  «  T:  e  T}.  Using  the  same  arguments  as  in  the  proof  of 
Proposition  1  it  is  easy  to  show  that 


tp(e,TP)  =  t®  -  iP(tP) 


From  Propositions  1  and  2 


c*(e)  =  1®  -  lP(H®) 
P 


where  1®  =  min{l*^(H)  :  H  s  ff®} 
P 


The  problem  (6)  is  a  relaxation  of  the  problem  (9),  which  implies  that 

1®  >  t®  and  now  from  (7)  and  (8)  we  have 
P  P 


'*NWWWVww¥\7»Tjvvrvnf  iT*  ir» 


c"‘(e)  >  tp  (e,TP)  +  iP(tP)  -  lP(H°)  =  tp(e,TP)  -  A(p) 


An  analogue  of  Lemma  1  may  be  also  proved  for  the  SHIP  and  the  S1TP  as  its 
relaxation. 

Some  comments  concerning  Lemma  1  are  necessary.  Two  special  cases  have 
to  be  considered: 


CASE  1°  -  when  there  exists  a  penalty  vector  p*  s  R  such  that 

A(p»)  =  0; 


CASE  2°  >  when  there  is  so  called  duality  gap  A  >  0,  where 

A  =  inf{A(p)  :  p  «  R*^}. 


In  the  Case  1°,  H°  =  arg  min{lP  (T)  :  ^  s  T}  and  from  Lemma  1  we  have 
the  following  inequalities  for  e  «  E: 


c'^(e)  >  t^,,(e,H°) 


c"(e)  >  tp,(e,H  ) 


In  section  3  we  will  show  that  bounds  for  c‘^(e),  c~(e)  provided  by  the 
inequalities  (10),  (11)  may  be  slightly  improved,  because  in  the  Case  1° 
stronger  inequalities  hold: 


c'*'(e)  >  t*(e,H°)  +  min{t\(u,H°) :  u  c  H°\{e}} 


c'(e)  >  t~,(e,H°)  +  min{t~M  (u,H°):  u  «  E\H°\{e}} 

p*  p" 
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In  order  to  use  inequalities  (10),  (11)  or  (12)  (13)  to  calculate  lower  bounds 
for  the  edge  tolerances  c^(e),  c~(e),  e  «  E,  in  the  Case  1°  two  problems 
have  to  be  solved: 


(i)  a  penalty  vector  p*  «  satisfying  A(p*)  =  0  must  be  found; 


(ii)  edge  tolerances  t'^M(e,H°),  t~«(e,H°),  e  «  E,  for  tne  SSTP  in  qP 

p  p* 

have  to  be  calculated. 


A  solution  of  problem  (ii)  is  described  in  section  3.  A  method  of  solving  the 
problem  (i)  is  discussed  in  section  4. 

In  the  Case  2°  bounds  for  c‘'’(e),  c~(e)  obtained  from  Lemma  1  are  weaker, 
because  A(p)  >  0  for  any  p  s  r'’.  Moreover,  Lemma  1  does  not  provide 
bounds  for  edges  belonging  to  (H°\T^)  u  (T^xH®).  In  section  3  we  will  prove 
a  theorem  which  specifies  bounds  for  c‘^(e),  c"(e)  in  this  case,  but  they  still 
may  be  weak.  Thus,  to  calculate  bounds  for  edge  tolerances  in  the  Case  2°  it 
is  required  to  find  a  penalty  vector  p,  for  which  A(p)  is  possibly  small  and 
the  cardinality  of  the  set  H°\tP  is  small  as  well.  This  problem  is  discussed 
in  section  4. 


3.  Edge  tolerances  for  shortest  spanning  tree  and  1-tree 

The  problem  of  calculating  edge  tolera.ices  of  the  shortest  spanning  tree 
has  been  addressed  in  several  papers  (see  e.g.  [1,  3,  12]).  In  this  section 
we  review  at  first  some  fundamental  facts  on  which  sensitivity  analysis  for 
the  shortest  spanning  tree  is  based.  Next  we  discuss  in  detail 
implementations  of  algorithms  for  finding  edge  tolerance  with  respect  to 
special  spanning  tree  which  is  also  a  Hamiltonian  path.  We  close  this  section 
by  proving  some  useful  result  concerning  relations  between  edge  tolerances  and 
lengths  of  spanning  trees. 


Let  T  be  the  shortest  spanning  tree  in  G  =  (V,E,C).  The  following 
well  known  proposition  (see  e.g.  [12])  formulates  necessary  and  sufficient 
optimality  conditions  for  T®: 


Proposition  3.  T®  is  the  shortest  spanning  tree  in  G  if  and  only  if  for  any 
e  «  E\T° 


c(e)  >  c(w)  for  w  «  U(e) 


where  U(e)  is  a  subset  of  edges  belonging  to  the  unique  path  in  T®  Joining  the 
ends  of  e.  D 


Denote  for  e  «  T 


W(e)  =  {w  «  E\T^  :  e  s  U(w)} 


The  following  fact  is  a  straightforward  consequence  of  Proposition  3: 


Proposition  4.  If  e  «  T^,  then  t~(e)  =  ®  and 


t'^(e)  =  min{c(w)  :  w  «  W(e)}  -  c(e) 


If  e  «  E\T  ,  then  t^(e)  =  ®  and 


t~(e)  =  c(e)  -  max{c(u)  :  u  «  U(e)} 


Proposition  4  provides  a  method  of  computing  edge  tolerances  with  respect 
to  T*^  by  finding  minimum  weight  edge  belonging  to  W(e)  for  any  e  «  T*^  and 


maximum  weight  edge  in  U(e)  for  any  e  ^  T  .  This  may  be  done  simultaneously 
by  using  an  auxiliary  graph  ([11])  called  transmuter.  A  transmuter  is  a 
directed  acyclic  graph  which  contains,  one  vertex  v(e')  of  in-degree  zero  for 
any  e'  «  one  vertex  v(e'')  of  out-degree  zero  for  any  e"  «  E\T*^  and 

arbitrary  number  of  additional  vertices.  Moreover,  in  a  transmuter  there 
exists  a  path  from  vertex  v(e')  to  vertex  v(e'')  if  and  only  if  e"  «  W(e'). 
It  was  shown  in  [10]  that  for  given  spanning  tree  T  in  G  a  transmuter 
containing  0(m  o(ra,n))  vertices  can  be  constructed  in  0(m  a(m,n))  time 
(where  a(m,n)  is  a  functional  inverse  of  Ackerman's  function  [10]).  Given  a 
transmuter,  a  labeling  procedure  was  described  in  [12]  to  compute  all  edge 
tolerances  in  0(m  a(m,n))  time  using  0(m)  space.  It  is  the  best  known 
complexity  of  algorithm  for  finding  all  edge  tolerances  with  respect  to 
general  shortest  spanning  tree,  although  there  is  some  doubt,  whether 
complicated  data  structures  used  may  lead  to  a  computationally  efficient 
procedure  (see  [10]). 

In  [3]  simpler  data  structures  were  proposed  to  compute  all  edge 
tolerances  for  (general)  shortest  spanning  tree  in  0(m  log  n)  time  using  0(m) 
space . 

In  [9]  two  methods  which  may  be  used  to  compute  edge  tolerances  were 
described:  the  first  has  time  and  space  complexity  O(n^),  the  second  has 

running  time  0(mn)  and  requires  0(m)  space. 

Any  of  the  methods  mentioned  above  may  be  used  to  calculate  edge 
tolerances  with  respect  to  the  shortest  spanning  tree  in  the  Case  2°  (see 
section  2).  But  in  the  Case  1°,  T®  is  a  particular  spanning  tree  which  is 
also  a  path  and  more  efficient  algorithms  may  be  proposed. 

Let  be  the  shortest  spanning  tree  in  G  =  (V,E,C).  Assume  that 

is  also  a  Hamiltonian  path  in  G  and,  moreover,  the  vertices  of  G  are  numbered 
in  such  a  way,  that  T®  =  [(1,2),  (2,3),...,  (n-l,n)}. 


Then  sets  U(e),  W(e),  e  «  E,  appearing  in  Proposition  4  are  defined  as 


follows: 


For  e  «  E\T^,  i.e.,  e  =  (k,l),  k=1,...,n-2,  l=k+2,...,n, 


U(e)  =  {(i,i+1)  :  k  <  i  <  1-1} 


For  e  «  T®,  i.e.,  e  =  (i,i+1),  i=1, _ ,n-1, 


W(e)  =  {(k,l)  :  1  <  k  <  i,  i+1  <  1  <  n,  (k,l)  i  (i,i+1)} 


Figure  1  illustrates  subsets  of  elements  of  edge  weight  matrix  C  for  ■  hich 
appropriate  minima  and  maxima  must  be  calculated  according  to  formulae  (16), 
(17)  and  (18),  (19). 


ss 


e  =  (i,i+l) 


If  the  graph  G  is  dense,  i.e.,  m  =  0(n  ),  then  the  following  simple 
labeling  algorithms  may  be  used  to  calculate  edge  tolerances 
t'^(e)  ,t"(e) ,  e  «  E,  in  OCn*^)  time  using  0(n  )  space.  Let  w(i,j)  «  R 

be  labels  defined  for  i  =  0,1, ...,n,  J  =  1,...,n,n+1. 

Algorithm  for  calculating  t'*'(e).  e  « 

Step  1  (Initialization)  for  i=1  ^n-1  ^  w{i,n+1):  =  <»; 
for  J:=2  ^  n  ^  w(0,J ):=“>; 

Step  2  (Labeling)  for  i:  =  1  ^  n-2  ^ 

for  J:=n  down to  i+2  ^  w(i,J):=  min{w( i-1 , j) ,  c(i,J),  w(i,j+1)}; 

Step  3  (Calculation  of  tolerances)  for  i:  =  1  ^  n-1  ^ 

t'*’(i,i+1)  :=  min{w(i-1  ,i+1),  w(i,i+2)}  -  c(i,i+1). 

Algorithm  for  calculating  t~(e).  e  «  E\T^ 

Step  1  (Initialization)  for  i:  =  1  ^  n-1  ^  w(  1, i+1 )  :=c( i, i+1 ) ; 

Step  2  (Labeling  and  calculation  of  tolerances) 


for  i:=2  ^  n-1  do 
for  j :=1  ^  n-i  do 

begin  w(J,J+i):=  min{w(  J,  J+i-1 ) ,  w(J+1,j+i)},* 
t~( j* J+i) •=«( J. Jf J+i) 
end 

If  the  graph  is  sparse,  then  more  efficient  algorithms  may  be  used  to 
calculate  edge  tolerances  with  respect  to  T®. 
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I  <r>t  4^1 


Consider  at  first  the  problem  of  calculating 


max{c(u)  :  u  «  U<k,l)} 


i* 
I' 


where  U(k,l)  is  given  by  (18).  In  order  to  solve  (20)  for  all  (k,l)  «  E\T^ 
efficiently  let  us  store  elements  of  set  U  =  {c(i,i+1),  i=1,...,n-1}  using  a 
data  structure  called  a  symmetric  heap  (see  [6]).  A  symmetric  heap  SH(U)  is  a 
directed  binary  tree  containing  one  vertex  for  any  element  of  the  set  U.  The 
vertex  v(i),  i=1,...,n-1,  of  SH(U)  has  a  label  c(i,i+1)  and  the  following 
properties  are  satisfied  for  any  k,l=1 , . . .  ,n-1 ,  k  1: 

If  c(k,k+1)  £  c(l,l+1),  then  there  is  a  path  in  SH(U)  from  the  vertex 
v(l)  to  the  vertex  v(k)  and,  moreover,  if  k  <  1,  then  v(k)  belongs  to  the  left 
subtree  of  v(l),  otherwise  v(k)  belongs  to  the  right  subtree  of  v(l). 

A  symmetric  heap  SH(U)  may  be  constructed  in  0(n)  steps  and,  as  it  was 
observed  in  [6],  any  particular  problem  (20)  for  given  k,l  is  equivalent  to 
calculating  the  nearest  common  ancestor  of  vertices  v(k),  v(l-1)  in  SH(U). 
But  this  problem  may  be  solved  in  0(1)  time  (see  [4])  if  a  preprocessing 
requiring  0(n)  time  has  been  performed.  This  means,  that  all  lower  tolerances 
t“(e),  e  «  E\T*^,  may  be  calculated  in  0(m)  time  and  0(m)  space. 

To  calculate  all  upper  tolerances  t'*’(e),  e  «  T®,  a  simple  algorithm 
requiring  a  sorting  of  values  c(e),  e  c  E\T®,  may  be  constructed  and  this 
problem  may  be  solved  in  0(n  log  m)  time  and  0(m)  space.  But  it  is  not  known, 
whether  there  is  of  linear  complexity  (0(m)  time  and  space)  algorithm  to 
calculate  upper  tolerances  of  edges  in  the  Case  1°. 

Edge  tolerances  with  respect  to  the  shortest  1-tree  T^  can  be  computed 
in  a  similar  way.  An  approach  is  based  on  simple  Proposition  5  which  is  an 
analog  of  Proposition  4  and  which  we  will  state  without  proof.  Let  be  an 


optimal  solution  of  the  S1TP  in  G  =  (V,E,C),  i.e.,  T  =  u  {(1,k),  (1,1)}, 
where  is  the  shortest  spanning  tree  in  the  graph  C) 

obtained  from  G  by  removing  the  vertex  1,  and  (1,k),  (1,1)  are  two  shortest 
edges  incident  to  the  vertex  1.  By  W^(e),  U^(e)  we  denote  subsets  of  edges  of 
Gj  defined  for  in  the  same  way  as  the  sets  W(e),  U(e)  for  T®. 

Proposition  5:  If  e  «  T^,  then  t~(e)  =  “  and 
t'^(e)  =  min{c(w):  w  «  W^(e)}  -  c(e). 

If  e  «  t'*’(e)  =  “  and  t~(e)  =  c(e)  -  max{c(u):  u  «  U^(e)}. 

If  e  «  E!j  =  E\E^\((1,k),  (1,1)},  then  t^(e)  =  >»  and 
t~(e)  =  c(e)  -  max{c(1,k),  c(1,l)}. 

Furthermore,  t"(1,k)  =  t"(1,l)  =  “  and  t'*'(1,k)  =  c(1,k)  -  min{c(e)  :  e  «  E^}, 
t'*'(1,l)  =  c(1,l)  -  min{c(e)  :  e  s  Elj}. 

D 

We  will  close  this  section  by  proving  a  result  which  establishes  a 
relation  between  the  edge  tolerances  with  respect  to  the  shortest  spanning 
tree  and  the  value  of  difference  between  the  weights  of  the  shortest  spanning 
tree  and  an  arbitrary  spanning  tree. 

Theorem  1 .  Let  T®  be  the  shortest  spanning  tree  in  G  and  T  be  an  arbitrary 
spanning  tree  in  G,  Then 

1(T)  -  1(T°)  >  max{  I  t'"(r),  I  t~(q)} 

r«T^\T  q«T\T*^ 


(21) 
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Proof.  Consider  two  subsets  of  edges:  R  =  T^\T  and  Q  =  T\T*^.  It  is  known 

(see  [5],  Theorem  1),  that  there  exists  a  bijection  from  R  into  Q,  such 

that  for  every  edge  r  «  R,  =  T®\{r}  u  {i|*(r)}  is  a  spanning  tree  in  G 

and  c(<)(r))  -  c(r)  >  0.  From  the  fact  that  is  a  spanning  tree  it 

follows,  that  i|)(r)  «  W(r)  and  from  (16)  we  have  the  inequality 

t'*’(r)  <  c(il»(r))  -  c(r)  and  further  1(T)  -  1(T*^)  =  ^  [c(\i)(r))  -c{r)] 

>  t^(r),  Similarily,  for  every  edge  q  s  Q,  T'^  u  {q}\{il)”  (q)}  is 

r«R  . 

also  a  spanning  tree  and  this  implies  that  iti~'{q)  «  U(q),  Now  from 

(17)  we  have  t“(q)  <  c(q)  -  c(4i"^(q))  and  finally  1(T)  -  l(T^)  = 

,-1. 


I  [c(q)  -  c(t|)"  (q))]  >  I  t"(q). 
qsQ  q«Q 


As  corollaries  of  Theorem  1  we  obtain  some  properties  of  edge  tolerances 
with  respect  to  the  shortest  Hamiltonian  path,  which  were  stated  without  proof 
in  section  2. 

Let  for  some  p  «  R*^,  H°  and  tP  be  optimal  solutions  of  the  SHPP  and  the 
SSTP  in  gP  =  (V,E,cP).  As  before,  tp(e,T^),  t~(e,T'^),  e  s  E,  are  edge 
tolerances  with  respect  to  T°  and  c^(e),  c~(e),  e  «  E,  are  edge  tolerances 
with  respect  to  H°, 

Theorem  2.  If  A(p)  =  0  and  =  tP,  then  for  e  «  E 


c'^(e)  >  f^(e,H°)  +  min  {t'^{u,H°)  :  u  «  H°\{e}} 
P  P 


.  ^wifc  b'l 
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Proof.  If  c'*’(e)  <  »,  then  according  to  (1)  we  have  c'*’(e)  =  lP(H®)  -  lP(H°), 

where  H®  =  arg  min  {1^(H)  :  H  s  H®).  It  is  easy  to  see  that  |H°\H®|  2. 

Obviously,  H°,  H®  «  T  and  now  because  H°  is  the  shortest  spanning  tree  in  gP 

and  e  «  H°\H®,  from  (21)  we  have  c'*^(e)  =  1^(H®)  -  1^(H°) 

>  t^(e,H°)  +  t^(u,H®)  for  some  u  «  H°\{e}.  The  proof  of  second  part  of 
P  P 

theorem  is  analogous.  0 

Theorem  3.  If  H°  |  T^,  then  for  e  «  H°\tP 


c‘^(e)  >  min{t“(q,TP)  :  q  e  EvT^Me}}  -  A(p) 
P 


and  for  e  «  T^\H° 


c“(e)  >  min{tp(r,TP)  :  r  s  T^Xle}}  -  A(p) 


Proof.  We  will  prove  only  (22),  because  a  proof  of  (23)  is  analogous. 
Consider  e  «  H°\T^.  If  c'^(e)  <  ”,  then  c'^(e)  =  1(H®)  -  1(H®)  and  there 
exists  a  spanning  tree  T®  which  is  the  second  shortest  spanning  tree  not 
containing  e.  Moreover,  1(T®)  <  1(H®)  and  c^(e)  >  1(T®)  -  l(T^) 

+  (1(tP)  -  1(H°))  =  1(T®)  -  1(tP)  -  A(p).  But  t|\tP  must  contain  some 
edge  q  «  E\TP\{e}.  Now  from  (21)  we  have  1(T®)  -  1(tP)  >  t~(q)  which 
implies  (22).  D 

Bounds  for  edge  tolerances  provided  by  Theorem  3  (and  Leauud  1)  may  be 
weak.  In  particular  cases,  values  of  right-hand  sides  of  inequalities  (22), 
(23)  and  (4),  (5)  may  even  be  negative,  which  means  that  trivial  bounds  are 
obtained.  Thus,  although  any  penalty  vector  may  be  used  to  calculate  c‘^(e). 


.'Tw  VJ 


c“(e),  it  is  desired  to  have  a  vector  p  which  gives  small  (if  possible  -  equal 
to  zero)  values  of  A(p)  and  |H°\tP(,  This  problem  is  discussed  in  the  next 
section. 


4.  Computing  of  penalties 


To  calculate  lower  bounds  of  edge  tolerances  with  respect  to  a 

penalty  vector  p  is  needed,  for  which  A{p)  and  |H°\Tp|  are  as  small  as 
possible.  If  the  duality  gap  A  is  equal  to  zero,  then  such  vector  may  be 
found  as  a  solution  of  equation  A(p)  =  0  and  this  guarantees  also  that 
|HO\tP|  -  0.  Otherwise,  one  may  try  to  solve  this  bicriteria  problem  by 


choosing  as  a  vector  p  such  a  feasible  solution  of  equation  A(p)  =  A,  for 
which  |H°\tP|  is  minimal. 

To  solve  A(p)  =0  two  attempts  may  be  considered: 

(i)  The  problem  rain{A(p)  :  p  «  r’’}  may  be  solved  exploiting 
properties  of  the  function  A(p)  (  A(p)  is  convex,  piece-wise 
linear  function  on  by  some  subgradient  type  procedure. 

(ii)  A  feasible  solution  of  A(p)  =  0  (if  exists)  may  be  calculated  by 
finding  a  solution  of  the  system  of  linear  inequalities  (24). 

The  later  approach  was  used  in  computer  implementation  and  it  will  be 
described  in  this  section. 
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Define  for  a  given  graph  G  =  {V,E,C),  P(C)  =  {p  «  r”: 


p(i)  +  P(J)  -  p(k)  -  p(k+1)  >  c{k,k+1)  -  c(i,j) 


for  (i,J)  «  E,  i  =  1,...,n-2 
j  "  i+2  ^ y  n 
k  =  i,...J-1} 


(24) 


Theorem  4.  Let  H°  =  {(1,2),  (2,3) , . . • , (n-1 ,n) }  be  an  optimal  solution  of  tne 
SHPP  in  gP  =  (V,E,CP).  Then  A(p)  =  0  if  and  only  if  p  s  P(C). 


Proof.  A(p)  =0  if  and  only  if  is  also  an  optimal  solution  of  the  SSTP 
in  gP,  i.e.,  if  necessary  and  sufficient  optimality  conditions  formulated  in 
Proposition  3  are  satisfied.  This  means  that  for  H°  inequalities  (14)  must 
hold  for  the  graph  gP.  But  for  the  spanning  tree  H°  the  sets  U(e),  W(e)  are 
given  by  (18)  and  (19),  and  now  it  is  easy  to  check,  that  if  the  inequalities 
(14)  are  formulated  for  H°  and  the  graph  gP,  then  we  obtain  a  system  of 
conditions  defining  P(C).  D 

The  number  S(G)  of  inequalities  defining  P(C)  is  of  order  0(mn).  If 
G  =  Kj^  (complete  graph  with  n  vertices),  then 


S(K^)  =  i  (2k-1)  [k(2k+1)-3] 


if  n  =  2k,  k=1 ,2, 
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Any  vector  p  belonging  to  P(C)  may  be  used  as  a  penalty  vector  to  compute 
lower  bounds  of  edge  tolerances  with  respect  to  although  different  vectors 
lead,  in  general,  to  different  values  of  these  bounds.  If  P(C)  =  ♦,  then  it 
means  that  there  is  a  positive  duality  gap  A. 

As  a  simple  consequence  of  Theorem  4  we  obtain  the  following  fact: 

Corollary  1 .  If  for  a  given  graph  G  there  is  zero  duality  gap  A,  then  the 
optimality  of  arbitrary  Hamiltonian  path  may  be  verified  in  a  polynominal 
time . 

Proof.  It  is  an  immediate  consequence  of  the  fact  that  P(C)  is  defined  by 
polynominal  number  of  inequalities  and  its  consistency  may  be  checked  in  a 
polynominal  time  by  linear  programming.  0 

Similar  facts  (which  we  will  give  without  proof)  hold  for  the  SHIP. 
Define  for  G  =  (V,E,C),  P(C)  =  {p  «  r": 

p(k}  -  p(2)  >  c(1,2)  -  c(1,k)  for  k=3,...,n-1,  (1,k)  s  E, 

p(k)  -  p(n)  >  c(1,n)  -  c(1,k)  for  k=3,...,n-1,  (k,n)  «  E, 

p(i)  +  p(J)  -  p(k)  -  p(k+1)  >  c(k,k+1)  -  c(i,j) 

for  (i,J)  «  E,  i=2, . . . ,n-2, 

J  ”  i "^2 ,...,n,  k=i,...,J“1} 


J.V  U’V  \.v  L'V’.Vi^’VVVffVVLV^.VVV-’ 
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Theorem  5 .  Let  R°  =  {(1,2),  (2,3) , . . . , (n-1 ,n) ,  (n,1)}  be  the  shortest 
Hamiltonian  tour  in  gP  =  (V,E,CP).  For  FP  to  be  the  shortest  1-tree  in  gP  it 
is  necessary  and  sufficient  that  p  «  P(C).  D 


5.  Implementation  of  the  method  and  conclusions 

The  method  of  calculating  lower  bounds  of  edge  tolerances  for  the  SHPP 
and  the  SHTP  described  in  previous  sections  was  implemented  for  IBM  PC  in 
Turbo  Pascal  3.0.  In  the  step  of  computing  of  edge  tolerances  for  spanning 
trees  and  1-trees  simple  O(n^)  labeling  procedures  mentioned  in  section  3  are 
used.  To  calculate  appropriate  penalties  an  approach  provided  by  Theorems  4 
and  5  is  used.  Penalities  are  computed  by  solving  linear  programming  problems 

T 

min{a  p:  p  «  P(C)} 


minla^p:  p  s  P(C)}  (25) 

Different  objective  vectors  a  «  may  be  chosen  and,  usually,  different 
penalties  as  well  as  different  lower  bounds  for  edge  tolerances  are 

obtained.  In  computational  experiments  a  =  (1, 1  )'^  or  a  =  (0,...,0)^  was 

mainly  used.  In  the  later  case  by  solving  (25)  an  existence  of  feasible 


solution  of  equation  4(p)  =  0  is  checked. 

To  solve  (25)  a  simple  specialized  version  of  the  revised  simplex 
algorithm  was  implemented.  As  problem  (25)  has  only  n  variables  and  large 
number  of  constraints  (for  example,  for  the  SHPP  in  n  =  40,  the  number  of 
constraints  exceeds  10000),  the  dual  problem  for  (25)  is  solved  and  column 


?iiiii>fiu  KtfiMTfi  nitmtii  fiawrrr  muTrrr  \tiiii 


generation  technique  is  used.  The  computational  experience  is  limited  to 
rather  small  sizes  of  problems.  In  Table  1  computation  times  in  seconds  for 
IBM  PC/XT  with  math-processor  are  reported.  These  times  do  not  include  input 
and  output  of  data.  All  test  problems  were  randomly  generated  as  planar 
Euclidean  SHTP. 

In  Table  1  n  denotes  the  number  of  vertices,  5  is  a  density  of  graph, 

T^.  ,  are  respectively  -  average,  minimal  and  maximal  times  of 

a’  min  max  ^ 

computing  penalties  (for  5  problems),  is  a  time  of  computing  edge 

tolerances . 

TABLE  1 


n 

5 

tP 

a 

xP. 

min 

xP 

max 

t 

X 

10 

1 

2.6 

0.5 

4.3 

0.1 

10 

0.3 

1.2 

0.3 

1.8 

0.1 

25 

1 

121.6 

58.8 

177.5 

0.7 

25 

0.3 

48.9 

40.6 

49.8 

0.7 

40 

1 

1572 

877 

2523 

1.9 

40 

0.3 

421 

316 

747 

1.9 

An  approach  described  in  section  2  may  be  used  with  different  relaxations 
of  original  problem. 

Let  (P):  min{f(x)  :  x  «  X}  denote  the  original  (primal)  programming 
problem  and  let  x°  be  its  optimal  solution.  Denote  by  (Rq)  a  relaxation  of 
(P)  parameterized  by  some  element  q  belonging  to  a  specified  set  Q: 

(Rq)  v(q)  =  rain{fq(x)  :  x  «  X^}. 


(For  example,  in  the  approach  described  in  this  paper  a  role  of  parameter  q  is 
played  by  the  penalty  vector  p  and  Q  =  R^). 

As  a  dual  problem  for  (P)  the  following  problem  may  be  considered: 

(D)  q»  =  arg  max  v(q) 
qcQ 

The  relaxation  (Rq»)  seems  to  be  a  good  candidate  to  provide  a 
sensitivity  analysis  for  x°  by  similar  approach  as  used  in  this  paper.  In 
order  to  apply  this  approach  one  must  be  able  to  answer  the  following  two 
auxiliary  questions: 

(i)  How  to  perform  a  sensitivity  analysis  for  the  problem  (Rq»)? 

(ii)  How  to  find  in  an  efficient  way  q*,  if  x°  is  given,  i.e.,  how  to 
solve  the  dual  problem  when  the  solution  of  primal  is  known? 

In  some  cases  an  answer  for  the  later  question  is  obtained  as  an  inexpensive 
by-product  of  solving  the  original  problem.  This  may  be  an  important  argument 
for  the  choice  of  relaxation,  because  as  numerical  results  reported  in  Table  1 
show,  in  this  approach  almost  all  computing  time  may  be  spent  on  solving 
problem  ( ii ) . 
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